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SUMMARY 


The flow of a compressihle fluid through a channel having locally 
supersonic regions is studied hy vising the Tricomi equation in the hodo- 
graph variables as an approximation in the sonic region to the equation 
of flow of an irrotational, inviscid gas. It is shown that this is 
equivalent to studying the flow of a gas having a pres sure- density rela- 
tion matching the Isentroplc relation to the third derivative at the 
sonic point. A one-parameter family of solutions of the Tricomi equa- 
tion is used which provides symmetrical accelerated-decelerated flows. 
The variation of this parameter alters the Mach number at the center of 
the throat, the velocity distribution and gradient along the center, 
streamline, as well as the shape of the channel, that is, the curvatvire 
of the bounding streamline. 

As specific examples, flows are computed having Mach number equal 
to unity and to 0.86 at the center of the throat section. Constant- 
velocity lines are plotted and it is found that the velocity gradient 
becomes zero at three places along each . streamline outside of a limiting 
streamline for values of the parameter greater than zero (M < 1 at 
center of throat section). For the p'arameter equal to zero (M = 1 at 
center of throat section), the velocity gradient along the streamlines 
and the curvature is discontinuous at all points of the two character- 
istics which meet the center streamline. 

Other solutions to the Tricomi equation are discussed which may be 
used to fomulate channel flows. The exact nature of these flows has 
not yet been investigated. 


INTRODUCTION 


In the investigations of the effect of free- stream Mach number on 
the lift and drag of an airfoil in a unlfoorm stream, it was found that 
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the lift did not begin to decrease or the drag increase until the free- 
stream velocity exceeded the critical free-stream velocity, that is, 
the velocity for which the local velocity of sound is just attained at 
some point of the body. Thus it appeared that, for free-stream velocities 
in between the critical free-stream velocity and the velocity for which 
the lift characteristics changed rapidly with Mach number, there might 
exist about the body a locally supersonic region which contained no shocks, 
although such shock-free flows have never been observed experimentally. 

A similar phenomenon occurs in the study of flows through a gradual 
constriction. When the exit pressure of a nozzle is decreased, the 
mean velocity at the narrowest cross section and also the total mass 
flow are increased. However, when the mean velocity at the constriction 
reaches the local velocity of sound, the mass flow reaches a maximum; 
and the flow then changes from an accelerated- decelerated flow to one 
in which the fluid is always accelerated in passing through the nozzle. 
Fiirther decreases in the exit pressure fall to Increase the mass flow 
but only introduce shocks into the supersonic side of the nozzle to 
satisfy the conditions of conservation of mass, momentum, and energy. 

One question which occurs regarding such a process is: Do there occur 

isolated regions of shock- free supersonic flow in the neighborhood of 
the walls before the mass flow reaches the TtinxlTnum and the flow becomes 
pvurely accelerated? 

This work was done at Brown University under the sponsorship and 
with the financial assistance of the National Advisory Committee for 
AeronaTitics . 


SYMBOLS 

?c,y Cartesian coordinates (z = x + ly) 

u,v velocity components in x- and y-dlrection, respectively 

7 ratio of specific heats 



Mach angle 

T = q^/2p (t = 1 corresponds to limiting speed) 
^ velocity vector 


<1 


magnitude of velocity vector 
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0 inclination of velocity vector vith x-axis 

= 2t(1 - t)"^ 


T2 


1 - (2p + 1 )t 
2t(1 - 


c local speed of sound normalized to speed of sound in 

still air 


p density of fluid normalized to stagnation density 

V specific volume 

p pressiire of fluid corresponding to density p 

Ps critical pressure (pressure at s = O) 

M Mach numher (q/c) 


s = 


2.1120s* 

1 + 1.230lte* 


s* = 


'( 213+1) -1 (1. 

2t 


dr 



qp velocity potential 

ijr stream function 

i|f« = oijjf ^ (1. • 4.5320 - 0. 84664s )i|r 

qp’ conj;igate potential to t' ("by equations (6)) 

$ symmetrical stream function given hy equation (19) 

stream function given hy equation (15) 
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EARLY lEVISTIGATIONS OF LOCALLY SUPERSONIC REGIONS 


The first attempt to construct a flow through a channel having 
locally supersonic regions was made hy Taylor in 1931 (reference l). 

He expanded the velocity potential in a double power series about some 
central point of a two-dimensional nozzle whose walls consisted of two 
arcs of circles. Taylor was able to obtain several flows which were 
symmetric in the velocity about the line of centers of the two circles. 
The values of the maximum velocity on the straight streamline for these 
flora were bounded by. a certain velocity for which the local velocity 
of sotind on the wall was J\ist barely exceeded. From this resiilt. 
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Taylor concluded that the potential flow through a channel with locally 
supersonic regions was impossible. Gortler (reference 2 ), however, 
felt that Taylor's conclusion was unfounded, for if more terms of the 
power series were taken and the boundary conditions on the wall satisfied 
to more terms in the series expansion, then the difficulty which Taylor 
encountered in determining the coefficients might not exist. Gortler, 
using much of Taylor's computations, satisfied one less relation of the 
boundary conditions than Taylor and thereby had two free parameters to 
vary instead of one. One of these parameters represented the velocity 
at the center axis of the throat of the nozzle, and the othdr, the 
velocity gradient along the center axis at this point. By varying the 
velocity gradient Gortler obtained a sequence of channels which exhibited 
locally supersonic regions near the channel walls at the constriction. 

As the velocity gradient was increased, the locally supersonic region 
gradually extended toward the axis. When the local velocity of sound 
was attained at the axis, the flow became identical to Taylor's unsym- 
metrical flow through a Laval nozzle. Gortler showed that the locally 
supersonic region cannot extend to the center axis if the velocity 
gradient is to be continuous along the center streamline. When the 
sonic velocity is just attained on the axis, but not exceeded, then the 
condition of continuity of the gradient will require a zero gradient at 
the sonic point. Under this condition the following theorem, proved by 
Gortler, shows that the sonic line must be a straight line and the 
curvatiore of the wall at the throat must be zero. 

Theorem ; Let an Inviscid flow of a compressible gas 
be given which is symmetrical about the y = 0 axis and is 
analytic in the neighborhood of a point x = Xq and y = 0. 

Further, suppose that the ^ocal velocity of sound is attained 
at this point of the center axis and that the velocity gra- 
dient vanishes; that is, u(xq,o) = 1 and ^^Xq,o) = 0. 

Then it follows necessarily for all y's in the nei^bor- 
hood that u(xQ,y) = 1, |^(^o^y) = ■’'(^o^y) = 

Flows having locally supersonic regions in which the sonic velocity is 
just reached on the center axis, therefore, will not be analytic in the 
neighborhood of the center-axis point of the throat. Such a flow has 
been computed by Tomotika and Tamada and a similar flow -having symmetry 
about both horizontal and vertical axes is given in this work. The 
velocity gradient in both flows has a jump at the center of the nozzle. 

The work of Tomotika and Tamada (reference 3)^ using an approximate 
gas law, and of Frank! (reference 4 ), Gortler (reference 2), Oswatitsch 
and Rothstein (reference 5)^ Fox and Southwell (reference 6), anfl 
others, using power series, already indicates strongly that solutions 
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■whicti have locally supersonic regions for the exact flow equations do 
exist. Recognition must he given, of course, tO' the lack of knowledge 
about the convergence of the power series, and to the fact that the 
Trlcoml equation which the Japanese authors use approximates the actual 
gas only in a small neighborhood of the sonic velocity. Yet, the fact 
that both methods give flows having locally supersonic regions clearly 
points toward a similar resiilt from the exact equation for an isentropic 
compressible, inviscid fluid. But the existence of these solutions by 
no means settles the matter. It is to be noted that "existence" is 
used in the mathematical sense. ActuEilly, all experimental resTilts 
(Liepmann (reference 7) and Ih'ankl (reference 4)) show that irrotational 
flow breaks down soon after the appearance of the lower critical Mach 
number (Tsien and Kuo (reference 8)). None of the experiments shows 
locally supersonic regions of any appreciable extent without shocks. 


GENERAL PROPERTIES OF A LOCALLY SUPERSONIC REGION 


Frankl (reference 9 ) and Nikolsky and Taganoff (reference 10) have 
made some general investigations of the properties of the locally super- 
sonic regions without shocks. Frankl (reference 9 ) showed that the 
symmetric flow over an obstacle which is symmetric about two perpen- 
dicular axes is \miquely determined by the boundary conditions on the 
portions of the arcs DFD’ and EGE' in figure 1 which lie outside of the 
region defined by the two pairs of characteristics BD and BE and B’D' 
and B*E* and the conditions at infinity. The segments of the body DE 
and D‘E’ between the characteristics then depend upon the rest of the 
body and the free-stream Mach number. Thus it appears that symmetric 
flow with supersonic regions does not exist for an arbitrary body, 
although for a given free-stream Mach number there might well be a 
large class of profiles for which such flows do exist. 

The results of Frankl can be shown more rigorously and easily by 
making use of the following uniqueness theorem proved by Tricomi 
(reference ll): 

In a domain bounded by a ctirve in the subsonic region 
with its ter mini on the sonic line in the holograph plane 
and two intersecting characteristics drawn from these ter- 
mini, the solution of the Tricomi equation is uniquely deter- 
mined by the values prescribed on the subsonic arc and on 
one of the characteristic segments boxinding the domain. 

The character of the flow through a nozzle having locally supersonic 
regions is shown in the holograph plane in figure 2. If the value of 
the stream function is given on the arc BCD and the segment of the 
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axis OAB and the characteristic ED, then the solution is given uniquely 
inside the region BCDECB. In the physical plane, this means that the 
solution is given uniquely in the region A'B'C'D'E’O’A’ (fig. 3 ) since 
the arc OA in figure 2 maps into the line 6 = 0 (O' A' in fig. 3 ) the 

line AB into the center streamline, and the characteristics DE and EO 
into the Mach lines D'E' and E'O'. The same thing can he shown for 
A'O'G'H'I' J'A' since symmetry is assumed. Since the values of 6, q, 
nnd T|r are known on the two Mach lines E'O' and O'G', then. the solution 
is determined inside the quadrilateral E'L'G'O'. Thxxs the flow field 
is entirely determined hy the boundary conditions on the arcs C'F' 
prid M'l', and by the velocity distribution on the center streamline JB' 
and on C'B' and I'J'. 

The fact that, for an Inviscld flow which has locally supersonic 
regions, the boundary cannot be chosen in an arbitrary manner actually 
has little significance in the flow of real gases. In such flows, there 
exists a boixndary layer near the walls and a pertinent question regarding 
locally supersonic regions may be stated as follows: Does the boundary 

layer adapt itself so that locally supersonic regions can occur without 
shocks, or are shocks always present in locally supersonic regions? 

Nikolsky and Taganoff have discovered some of the properties which 
a locally supersonic region in an inviscid flow free of shock waves 
must have. Briefly these properties are as follows: 

(1) The slope of the contour at any given point in the locally 
supersonic region is the arithmetic mean of the slopes of the velocity 
vector at the points on the transition line which are the termini of 
the characteristics originating from the given point of the contour. 

(2) If one moves along the transition line in such a way that the 
region of subsonic velocities lies on the left, then the velocity vector 
will always rotate clockwise. 

(3) If; in the supersonic zone, a segment of a characteristic of 
one family is given such that the characteristics of the other family 
originating from points of this segment extend to the transition line, 
then the Inclination of the velocity vector and its magnitude are 
monotonlc functions along the given segment of the characteristic. 

(k) The magnitude of the velocity along a rectilinear section of 
the profile in a supersonic zone must decrease in the direction of flow. 

( 5 ) Prandtl-Meyer flow between two characteristics of one family, 
with the characteristics of the other family rectilinear, cannot be 
completely realized up to the treinsitlon line in the finite region. 
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(6) CliEiracteristlcs of the flxat family which originate from points 
of a rectilinear section of the profile cannot extend to the transition 
line^ hut must instead terminate in a compression shock’. 

Nikolsky and Taganoff also discuss the possihillty of the Mach lines 
for min g an envelope in the supersonic zone. They show that if the solu- 
tion to the flow over a hody is found which satisfies the hoimdary con- 
ditions then the Jacobian of the transformation from the holograph plane 
to the physical plane cannot vanish on the boundary, since the curvature 
of the streamline must be zero when J = 0. By using the property of 
statement (3) above they show that if the Jacobian does not vanish on 
the boxmdary then it cannot vanish in the locally supersonic region in 
the flow field outside the body. They give as the criterion of break- 
down the relation 


— ^ = q tan an 
d0 


where q is the velocity magnitude, 9 the angle of inclination of the 
contour, and the Mach angle associated with the local Mach number. 

What Nikolsky and Taganoff failed to notice, hoirever, is that the above 
equation is the condition that the streamline representing the body be 
tangent to a characteristic in the holograph plane. If the derivatives 
with respect to 6 and q of the stream function are finite, then the 
streamline will have infinite cxirvatiire at the point where the relation 
above is satisfied. What Nikolsky and Taganoff have actually shoim may 
be stated more precisely in the foUpwing form; 

Theorem ; If, along the portion of a streamline lying in 
a local supersonic region of flow, there exists no point at 
which the curvature becomes infinite, then there cannot exist 
any limit line inside the supersonic zone lying between the 
given streamline and that streamline which is jiist tangent 
to the transition line. 

However, if one considers the continiiation of the flow pattern then one 
may find a streamline for which the curvature is Infinite and hence the 
Jacobian J = 0. Nikolsky and Taganoff were unable to show that J = 0 
on the transition line was an impossibility under the conditions of the 
theorem above. Friedrichs (reference 12), however, shoired that, imder 
the hypotheses of the theorem above, infinite acceleration and infinite 
curvature of the streamlines could not occtir on the transition line 
between the two streamlines mentioned in the theorem. In investigating 
the proximity to breakdovm of a potential flow over an obstacle, it 
seems advisable to study the continuation of the flow into the region 
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inside tiie body and investigate how near to the body streamline lie the 
cusps of a limit line. If limit lines lie near the body, then it is 
probable that a slight change in the boundary 'or an increase in the 
free stream would result in an inability to satisfy the boundary 
conditions . 


APPROXIMATE HODOGRAPH METHOD FOR TRANSONIC FLOW 
INVOLVING SOLUTIONS TO THE TRICOMI EQUATION 


In this discussion an approximation to the compressible-flow hodo- 
graph equations will be used in which the Tricomi equation replaces the 
usual second-order linear equation for the stream function. In the 
notation of reference 13, the equations 




( 1 ) 


where T]_ = 2t(1 - t)”P and T2 = 


1 - (2p + 1 )t 1 - m2 


represent an 


2t(1 - ^ 

exact linearization of the equations of irrotational motion of an 
Inviscid, compressible fluid. If a new vsLriable 


I* = 


‘(2P+1) 


-1 


(i - t)P 


2t 


dx 


UT. 


is Introduced into equation (l)., the eqiiatlons 




= -^3* 


> 




(2) 
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are obtained, 
equation for 


wbere 03*( s*) = 
■>lr results when 



Then the following differential 


q) is eliminated from equation (2); 


ilTg^g* + aft(s*)tee = 0 


(3) 


The general forms for two simultaneous differential equations in 
two independent variables which reduce to the Tricoml equation have been 
developed by Loewner (reference l4). One such form will be considered 
here. Let ijr and q) be defined by the relations 


^ -Ps = 


W 


and let a transformation of the dependent variables be given by 


i}f' = 00]_ilr 


<P's = --^<Ps 


1 

<P'e = - + HT ♦ 


(5) 


On eliminating q)' from eqxiation (5)^ i't can be seen that, for com- 
patibility with equation (i^), the quantity 0 ^ must take the form 
(Op = as + b. i'fhen equations (5) are substituted into equation (4) the 
results are 
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Elimination of cp' from the above equations leads to the Tricomi 
equation 


"1^’ss + ^^'ee = 0 


( 7 ) 


Now equation (4) may be written in the form of equation (3) by the 
introduction of the independent variable s* defined by 


s* = 


ds 




+ Constant = 


b(as + b) 


if 


off- = 


b^s* 


(1 - abs*)5 


The function ts* can be made to approximate 


m2 


in the neighbor- 


hood of the velocity of sound best by determining the two constants so 


that the first two derivatives of oj* and 


m2 


coincide. From 


this, with y - 1.^, the res\ilts 

4 


00* = 


9.teo4s* 


(1 + 1.2304s*)5 


s = 


2.1120s* 

(1 + 1.2304s*) 


(8) 


and 


i = (1.45320 - 0.84664s )"H’ 


( 9 ) 


are obtained. The relation between b 
and tabulations of s, q, l/q, l/p, 
tables I and II. 


and s* is sho-vm in figure 4 
s*, and 00||^ are given in . 
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EXACT EQUATION OF STATE LEADING TO TEE TRICOMI EQUATION 


This reduction^ which consisted essentially of changes in the 
dependent and Independent variahles and an approximation of 

to the second-order term in its power-series expansion in terms of s*^ 
may he regarded as producing the exact equation of flow for a gas having 
a somewhat altered pressure-density relation. ^This would correspond 
to viewing the Karman-Tslen or the Chaplygin approximation for the holo- 
graph relations as the exact equations for a gas having a pressure- 

density relation of the form p = A + Using Bernoulli’s equation in 

differential form, the equation 



( 10 ) 


is obtained, where v = l/p. If the approximation oj’^(s-*-) for — ^ — 

is now \ised, a differential equation for v is obtained whose solution 
yields the relation for p as a function of s* for the new gas. 
Although the equation is nonlinear, it is of the Riccati type, and can 
be solved either by expressing v as a power series in s* or in closed 
form by making a suitable substitution for the dependent variable. Since 
the behavior of the flow near s* = 0 is desired, the power-series 
expansion was developed and it was foimd that the specific volume could 
be expressed as 


i = 1 . 577 ^ - 2.1<883s-)t + 8.635 Ms*)^ - 30.175(s*)3 + 88.024(s*)^ - 

223.32(s*)5 + . . . (11) 

The first four terms agree with a similar expansion of l/p for the 
isentroplc pressure-density relation with 7 = 1 . 4 . 

Using the relation between s* and q (equation (U)) and a 
Maclaurin expansion for l/q, the series expression for the velocity is 

i = 1.0954 + 1. 72808-’^ i.7199(s*)2 _ 3.93^(b*)^ 6.7097(s*)5 - 

9.5639(s*)^ + . . . 


( 12 ) 
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From this and Bernoulli '-s equation an expression may he obtained for the 
pressure as follows: 


P - Ps = 0.8333s* - 1.31^5(s*)2 + 0.6912(s*)3 _ 0.38l6(s*)^ + 

1.0939(s*)5 - i.39^(s*)6 + . . . ■ (13) 

Equations ( 11 ) and (13) then constitute the pressure-density rela- 
tion leading to the Tricomi equation. These expressions give agreement 
with the \isual isentropic relations up to the third derivative of p 
with respect to p at the sonic velocity. 

However, the use of a power series to obtain l/p is hampered by 
the fact that the series converges very slowly. As an alternative 
method, the Riccati equation may be solved in closed form. As usual 

this involves a substitution of the form v = i = — in equation (lO), 

P u 

which yields 


(l + cs*)5u" = d X 8*U 


where c = 1.2304 and d = 9*4204 are substituted for convenience. 
Letting 


u = (l +’ cs*)w(ti) 


where 





1 + cs* 


leads to the differential equation w^^ - qw' = 0. This is the well- 
known Airy equation and a linear combination of its solutions together 
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with houndary conditions on l/p, that is, at s* = 0 and i = l.^TJhh, 
provides the result 


1 ^ 1.230^ 

P 1.+ 1.2304s* 


where 


(.1 * l.a304s.)2 . 0.1963t^/3i_^^^(j^ 

I = / 1- 6118s* 

\i + i.23o4s*y 


and Iy(5) is a Bessel function with imaginary argument in the notation 

of Watson (reference 15 )♦ Hie quantity becomes imaginary for super- 
sonic values of speed (s* < O), and if 'g is defined by 

5 = ( — the relation for l/p in the supersonic range is 
\1 + 1.2304s*y 


1 ^ 1.2304 

P 1 + 1.2304s* 


2.4176 

(1 + 1.2304s*)^ 


~ 1/3 ^ ^1/3 ~ 

“I J1/3U) - 0.19631 J-i/3(0 


where Jy(5) is the Bessel function of the first kind. The func- 
tions q and M may also be conqnzted in the same manner and they sire 
found to have the form 
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t 


<1 = 


-0.1667 


(1 + 1.2304s*)[|^/\/3(|) - o.19635^^\_i/3(0 


= 1 9-4204s*p^ 


(1 + 1.2304s*)^, 


Thus again the state qiiantities have been found for the present 
hypothetical gas In tenns of s*. This pressxire- density relation has 
been compared graphically with the usual Isentroplc relation for 
7 = 1.4 (see fig. 5). The curves are almost coincident in the range 
P = 0.42 to p = 0.72. This corresponds to a Mach n^mlber range of 1.35 
to 0.85, For values of s* near zero, the ciorves plotted for Mach nvnn- 
ber as a function of s* for the two pressure-density relations agree 
quite closely (see fig. 6). For large values of s*, however, the 

approximation curve departs from the usual isentroplc curve p = p^/7, 
since the slope of the hypothetical pressiire-velocity relation decreases 
more rapidly t h a n that for the isentroplc curve (see fig, 5). 


HOPOGRAFS SOLOnONS YIELDING ACCELERATED- 
DECELERATED CHANNEL FLOW 


Solutions for the Tricorn! eqmtion (7) which produce an unsymmet- 
rlcal accelerated-decelerated channel have been given by Tomotika and 
Tamada (reference 3), Falkovlch (reference 16), and Ehlers (refer- 
ence 13). Consider a solution which has the form 


to(s;0) 







(14) 


Tomotika and Tamada noticed that this function remained a solution 
when (0 - i\) is substituted for 6, In fact, if the function 


)^l(s,e) = Re 


w 

pe - ^ - (e - 




e - ±x)^ + i s3 + (0 





(15) 
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is considered, a stream function results which has a hranch-point sin- 
gularity of the one-half power at ® 

branch point in the subsonic stream may he said to "generate" a channel 
flow with locally supersonic regions having a Mach number at the nar- 
rowest cross section determined hy X. Of the other properties of the 
function it may he observed that,, when s > Sq and 0 = 0, then 

= 0. It would appear that ilr^Cs,©) has a line of one-third-order 
branch points along s = 0. However, this branch point may be elimi- 
nated by writing the function in the form 


Jl/3 ^ 


and 0=0. This 
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where 



Since does not vanish (the choice of sign for the one-half -order 

branch point is such that this is insured) it can be seen that the sin-' 
gularities along s = 0 are removable. The observation may be made 
that 

tl(s,0) + ti(s,-0) = $(s,0) (19) 

is a symmetrical function of 6, and it is this function which shall be 
used as the stream function. The function f has been plotted in fig- 
ure 7 for X = 0 . 050912 . This corresponds to s = 0 .I 8 at the center 
of the channel or a maximum Mach . number on the center streamline of 0.86. 

For sufficiently large negative values' of s ( supersonic flow ) , ^ 

becomes zero at a point above the 0=0 axis as well as at 0=0. 

This Indicates the flow in the supersonic region will have streamlines 
along which the velocity gradient changes sign more than once. This 
can be shown from the following analysis. A series expansion of Ifg 
for constant s about 0 = 0 '.has the form . , , - 

■ (We=o ® ■ ■ 
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From the hodograph, it Is apparent that for s > 0.l8 the derivative 
$0 Is positive for 6 > 0. Along a streamline for which does 

not vanish for a value of 0 > 0 , the value of ^qq then will he 
positive in the neighborhood of the line 0=0 since ^ 0 for 
0^0. The limiting streamline of those along which the velocity gra- 
dient does not change sign more than once is the streamline which passes 
through the point on the s-axls for which vanishes. This tras 

found to he f = 0.762 for s = -0.132. For s < -0.132 and hence 
$ > 0.762 the value of for 0 > 0 in the neighborhood of the 

line 0 = 0 is negative; for s > -0.132 and $ < 0."{62, is 

positive. (See fig. 8 (a).) The resulting effects on the constant- 
velocity curves in the physical plane are sho™ in figure 8 (h). It is 

apparent from the figure that ^ = 0 does not occur along the charac- 

teristics passing through s = 0 and 0=0. This remark is made at 
this point since in the limiting case as X, — 0 this characteristic 
does become the locus of the changes in sign of velocity gradient. This 
is discussed in more detail in the section Limiting Case 'VJhen Sonic 
Velocity Is Just Reached on Center of Streamline. 

The breakdown of isentropic, potential flows is usually associated 
with the appeanance of limit lines. Some remarks have already been made 
about cr iter ions for breakdown, b^xt here it will be sufficient to remark 
that ToUmein (reference 17) and Tsien (reference 18) have shora that 
Inviscid potential flows cannot be continued across limit lines. Limit 
lines occur when the Jacobian of the transformation from the holograph 
to the physical plane 



( 20 ) 


becomes zero. This may only occur for s < 0 since s, and 

are real and the Jacobian is always positive for s > 0. Geometrically, 
a necessary and sirfficient condition that limit lines occur is that the 
streamlines and the characteristics be tangent in the holograph plane, 
provided the derivatives fg and $s remain finite. As a matter of 
fact in this case the limit lines are the Ioctis of the points of tangency. 


The holograph of the function $ was tested graphically for limit 
lines by seeking points of tangency of the streamlines using a "portable" 
characteristic curve, drawn on tracing paper and moved across the super- 
sonic field. From figure 7 it is apparent that no limit lines occur for 
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that portion of- the b ,6 , field which has Taeen computed. Using this 
test, it is not certain that limit lines do not appear somewhere else 
in the supersonic portion of the flow, hut since, the choice of hounding 
streamline has so far heen left open this is not an important prohlem. 


CALCULATION OF PEfSIGAL COORDHtATIB 


To obtain the coordinates of the flow in the physical plane the 
condition of irrotationality eind continuity may he put into the 
following form: 


dq) = u dx + V dy 


— dijr = (-V dx + u dy)i 
P 

These equations may then he written as 



Expressing cp and ijr as functions of s and 0, and then replacing 
the partial derivatives of qp with those of ijr hy equation (4), the 
following differential expression for the complex position coordinate 
is obtained: 


dz = + i te) de + ^so>L% + | tg) ds (22) 

With this formula it is possible to show that if is an even . 
f inaction of 0 the flow pattern in the physical plane is symmetric. 

For this purpose, it is first of all necessary to prove that the line 
0=0 in the holograph plane maps into a straight line in the physical 
plane. Setting 0=0 and d0 = 0 in equation (22) leads to 



ds 


(23) 
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The stream function $ in equation (i9) is an even function of 6 } hence 
^ = 0 for 9=0. Then 

dz = — ds = — d^ (24) 

pq^tl6=o Je=o . 

This is p\arely imaginaj^yj and 9 = 0 then is seen to map into the 
y-axis. If the integration for the streamlines is carried out along 
each streamline from the point 6=0, the obvious differential relations 


dilr = 0 = i)rg ds + Tjrg d9 


ds 



lead to the simple formula, where ijr- = Constant, 


z = 



(25) 


Since (see preceding section) the Jacobian of the transformation from 
the s,9- to the x,y-plane is never zero, and since it may be represented 
as 



T}f=Constant 
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or 


X = - 


pq cos 6 J 


d0 


ilf=Gonstant 


y = - 


Pe 


uo 


pq sin 9 J 


d0 


■»|r=Coiistant 


Now for f in equation (19)^ the relation p(0)q(0) = p(-0)q(-0) along 
any single streamline, since fg and $ are even in 0, implies that the 
same value of s occurs for 0 and -0. Thus for x .the integrand is 
a product of four even functions of 0: J, pq, fg, and _ cos 0. There- 

fore, X Itself is odd. In a similar' manner, y is sedn to be an even 
function of 0 . Since x is an odd function of 0 for a given value 
of s and y is an even function of 0, then the channel is symmetric 
about the line 0=0. 


The mapping from the holograph plane back to the physical plane 
may be^accon^jlished by integrating equations (21) or (22) along either 
constant ijr or constant s lines. Since either one of these inte- 
grations must be carried out by numerical or graphical methods, the 
choice depends on the method involving the least amount of labor,' 

If the integration is performed along a streamline, then from equa- 
tion (2l) the physical coordinates are found by integrating 


dz = i e^® dcp 
1 


If cp were known as a function of s and 0 along each streamline 
this would yield a very simple method for computing the flow. However, 
the differential dcp may be replaced by the differential of cp' with 
an additional term involving ijc and Multiplying the second of 

equations (5) by ds • and the third by d0 and adding the two relations 
lead to 

T - dOi 

dcp' = - — dcp + — - Tjr d0 
ds 
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Thus the physical coordinates are found hy integrating 


dz = - 


^ 

I 


+ e^®ilr d0 


(26) 


along ilf = Constant. 

Now the function cp* can he computed from equations (6) when ijr' 
is known. Differentiating with respect to s gives 


Tjr ' 
^Is 


= Re W 


i j(0 - il)2 + ^ s3 - (e - il^ 


-2/3 1 


2s^ 


3^(0 - iX)2 + ^ s3 


(0 - 


iX)^ + ^ s^ + (0 - IX) 


-2/3 


2 s^ 


3^(e - + I s 3 


This expression may he simplified by multiplying the first term inside 
the braces by \mity in the form 


[( 


(0 - iX)^. + ^ + (0 - iX^ 

(0 - 1X)2 + i s3 + (0 _ iX) 


2/3 


2/3 


and by multiplying the second term by a similar factor with the posi- 
tive sign before (0 - IX) changed to minus. This leads to 
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Since ^is' ~ 90'^ potential cp’ 'becomes 



The arhitraxy f\mction of s results from the integration with respect 
to 0. That f(s) =0 will 'be readily apparent for the functions ’ 
and qp' must'ohey the relation 


9s' = ~^e 



2h 


MCA TN 2547 


Thus, differentiating equation ( 28 ) with respect to s yields 



Employing the same procedure as on ilfg' leads to 



The expression on the right side cau be readily recognized as -s\}rg’ as 
required. The potential, cp’ corresponding to \lr’ given by equation (17) 
is then 



/ 
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If for constant s the function is an even function of 6, 
then cp' is an odd function. This is apparent from equation (6) above, 
since differentiation of an odd function of 0 with respect to 0 
resialts in an even function. Th\is, the conjugate potential corre- 
sponding to f is . • 


= cp'(s,0) - qp'(s,-0) 



3 







(30) 


Here use has been made of the fact that is an even fimction of 0. 

In order to perform the actual integration, the streamline and con- 
stant qp* lines must be plotted in the hodograph plane. Then for each 
streamline the q^aantities (axi/q) cos 0 and (coi/q) sin 0 are plotted 

against and also against 0. The integrations then are performed 
with the aid of a planimeter. 

A symmetric channel flow has been computed for Sq = 0.l8 
(X = 0 . 0509117 ), which corresponds to a Mach number of 0.86 at the 
center of the throat section (see figs. 8(a) pru^ 8(b)). 


VELOCITY DISTRIBUTION ON CENTER STREAMLINE f = 0 


The streamline $ = 0 must be treated separately. The position 
coordinates are computed by integrating equation (22), which for this 
case reduces simply to 


cbc . 
ds 


SOI, 2 
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since $g ~ ® nnfl d.0 = 0. iEtie derivative with respect to 6 is foimd 
hy differentiating the stream function in the form of equation (15). 

This leads to 


Re 

tt 


-n 

1/3 

J (0 - 1X)2 + 1 s3 _ (0 - iX) 


\ 

|(0 - iX)2 + ^ s3 


^(0 - iX)2 I (0 - iX) 

' ■ >> 

^(0 - il)2 + I s3 

- iX)^ +-^ s3 - (0 - 
^(0 - iX)^ + ^ s3 



(31) 


j(0 - 


IX)' 


ff 


(S)*"' 


4 s3|J(0 - ix)2 




(0 - IX^ 




(32) 


Finally, 


^10 = - 


1 

30 ^ 0 , 


cos 



(33) 
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The quantities a, cu, and ft are defined in equation (18). If 

0 goes to zero through positive values with s > then 


CD = 2jt and 


“ 3 5^1 (f - 3 ) 


Where 5 = tan“^ ^ and Oq = s3 - X^. If 0 goes to zero through 

/q p\V3 

negative values of 0 for ® ^ ^ 


then CD = 0 and 


tisf-o) = - 35^(5 cos I 


The derivative of the stream function for the symmetrical flow then 
becomes 

J0 =^t(e) + ^t(-0) 


80 


= te(e) - '^gi-e) 


which for 0=0 reduces to 


^0 = 


3cdicXo\ 


t- -’rE" 8 - 1 ) 


6 

+ cos — 

3 




. J|_ ^1/2 


(Ditto 


B ' cos 


(i-l) 


(34) 
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Witti this expression substituted for fg, the physical coordinates 
are found by integrating 


dx 


' ^( 1 ) 


Va 

qoto 


cos 


(f-l) 


ds 


( 35 ) 


Unfortvinately, = 0 at s 


= (i^) 


2/3 


and the quantity on the right- 


hand side becomes infinite so that numerical integration in this form is 

impossible. This can be remedied by noting that ^ ^ With 

ds 2 soq 

this substitution for Oq, equation (35) yields 


dx 

d8 




(36) 


The integration can be readily perfomaed by plotting the ri^t-hand 
side against S- and -using a planimeter. The velocity distribution 
along f = 0 as -vrell as along the other streamlines is plotted in fig- 
ure 9 for X = 0.050912. 


CDRVIS OF CONSTANT VELOCITY 


The curves q = Constant in the physical plane are found by inte- 
grating the relation 


dz 

d0 



q pq 



q=Constant 


( 37 ) 


To study the nature of the curve in the neighborhood of the vertical 
axis of symmetry x = 0 (0=0 and s < O.I8) it is convenient to use 

a po-wer-series expansion about 0 = 0. The coefficients of such a series 
are readily obtained by differentiating successively the expression on 
the right "vrith respect to 0 and then setting 0 = 0. The expansion 
for z becomes 
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since and fgg are zero at 0=0. The parametric representation 

of q = Constant in the neighborhood of 0=0 then has the form 


X = A0 + B03 + 0(0^) 
y = C0^ + 0(0^) 

The slope of the line q = Constant at the axis of symmetry is 

^ 2 G 0 + 0(03) 

^ ^ A + 0(02) 

d0 

nnd ^ = 0 at 0 = 0. Now the curvature is given hy 
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wMch for 0 = 0 reduces to 

a 2 \po^^s^^=0 

The sign of the curvature at the point 0 = 0 depends upon the sign 
of C which in turn depends on the sign of $gg. Now, the second 

derivative becomes zero for s = -0.132 and 0 = 0. For 

s > - 0.132, the quantity ^qq is greater than zero and the curvature 
of the s = Constant lines is convex to the centeb streamline of 
the channel (see fig. 8(b)). For s < -0.132, the quantity is 

less than zero for 0=0 and the curvature of the constant s line 
is concave to the center streamline of the channel. . The curves 
s = Constant for O.I8 > s > -0.132 show a minimum along x = 0, 
while the curves s = Constant for s < -0.132 exhibit a maximum at 
X = 0, and two minlmums symmetrically placed about the y-axls. Along 
the streamlines which cut these latter curves, the fluid is accelerated 
to a TtiflvlTTmm velocity and then is decelerated slightly to a minimum 
velocity at the throat. The fluid is accelerated to the maxinnim velocity 
before being decelerated finally to subsonic velocities again. This 
behavior is apparent from the velocity distribution along the streamlines 
$ = 0.832 and 0.920 in figure 9. 

The cxirves s = Constant for s >0.l8 which- Intersect the x-axls 
at 90*^ are concave to the incoming flow in the neighborhood of the 
straight streamline at the entrance side of the vertical axis of sym- 
metry and convex to the incoming floi^ on the exit side. Further from 
the center streamline and toward the wall, however, the cxurvature changes 
and the constant-velocity lines become convex to the incoming flow on 
the entiance side and concave to the incoming flow on the downstream 
side of the vertical axis of symmetry (see fig, 8(b)). 

Since the stream function has a branch point at s = 0.I8 and 
0=0, the above anadysis cannot be applied to the curve of constant 
velocity s = O.18. To investigate the nature of this curve in the 
neighborhood of the center streamline, it is necessary to obtain the 
expansion of eqxiation ( 37 ) for small values of 0. The derivative 
is found from equation ( 33 )- Using the fact that fi(-0) = jt - U(0) 
and cd(-0) = 2« - cu(0) leads to 
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*6 = 


JL 


cos 


(R - . s\ 

U 2 6 ) 


© 


1/3 


Pi 


cos 


fo. + 

V3 2 6J 


(39) 


where il, a, and cd are defined in equations (l8). In a similar way 
by using equation (27) for il^n_’(e) the following equations are 
obtained; 


Js 


$ = 2_ ^ 

® 9cJi a, 


gfV (f - I-I) 


Pl^ 






- —f 

o>i 


(ho) 


Letting s — 0.l8 in equations (39) and (^) leads to some simplifica- 
tion since for this case ^ Accordingly, for values of 0 

positive and neeir zero, 


a = 6^/2(e2 + (2 x 0)V2 ^ o(e5/2) 


. 


CD = tan" 


" k * 


+ X 


fl = tan" 




. g,^ + 

’ll 2 


02 


- 0 


Pi = 




+ 6‘ 


r-flap 


= xV3 


5 ^ iC - m * 


(^ 1 ) 
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The approximate values of the above quantities lead finally to 



= 0.729020" 1/2 1; + 2.7280 -1-0(0 3/2) + . . ] 


m 


ts 


^x 1/6/4\V3 g_i/2 

2 [9) (02_ 


2 a 




1 

o^X 


0 l/2 


= -0.30930 


- 1/2 


5.8770^/^ -I- 0(03/^) 


(i^ 3 ) 


With these expansions and equation ( 37 )^ the slope of the curve s = O.I8 
in the neighborhood of 0=0 is given by 


^ = 1 . 9 l ^5 - 39.530 + 0(02) 

dx 


The numerical value of the slope at 0 = 0 is finally 


i - 1 - 9 “ 


Nou the curvature is proportional to the second derivative 
for the present case is given by 


d^y 

^ which 
dx^ 


d0 \^J dx 


0^/^ -1- 0(e) 


The curvature of the s = O.18 line is seen to be zero where it cuts 
the center streamline. 



MCA TN 25l^7 


33 


LIMITING CASE WEEN SONIC VELOCITY IS JUST REACHED 
ON CENTER OF STREAMLINE 


Setting X = 0 in the stream function in equation (15) leads 
directly to the stream- function considered for the flow through a Laval 
nozzle in reference 13. From this solution it is possible, however, to 
construct a channel having locaJLly supersonic regions hy the appropriate 
choice of branches. Tonratika and Tamada (reference 3) constructed a 
channel having M = 1 at the center streamline using a solution in 
the ijr,q)-plane. Their solution for this limiting case can be shown to 
be the same as eqmtion (l^i) with X = 0 when expressed in the hodo- 
graph variables. Now equation (6) when expressed in cp' and ijr’ as 
the independent variables becomes 


m 


If s and cp’ were replaced by -s and -cp' the equation takes the 
form given by Tomotika and Tamada (reference 3) for K = i. Their 
solutions for equation (44) are then 




St’ = -®®cp’ 


I = -z - 2(t*)^ 

I = 2t' jz - 2cp’ + 


(45) 


where z must satisfy 


(z - 2 pi)^(z + |i) = a3 


and n = (t')^ - cp' . When a = Q it is quite simple to eliminate z 
and cp' from the equations above and to obtain t’ as a function 



of s snfl 0. Setting the factor z - 2ij. equal to zero and elimi- 
nating cp’ and z with the aid of equation (45) leads to 


ir)- 


3 - lift 

i6 


s + 


3 0 


= 0 


The resulting function ''I'’ is 



1 



12 t i s3)^/3 

'I • 9 / 


Similarly, the relation z + p. = 0 yields 


( 46 ) 


(t‘) 


3 



36 

8 


0 


whose solution leads to the stream function 



Equations ( 46 ) and ( 47 ) express the two branches of the stream function 
which go to maike up the accelerated-decelerated flow thro\igh a nozzle. 
Equation ( 46 ) may he obtained from equation ( 47 ) by inserting unity in 
the form of e^^^ inside of the parentheses on the right-hand side of 
equation ( 47 ) and taking the real part. Eor the flow computed by 
Tomotlka and Tamada, figure 10 shows the regions of the holograph plane 
in which the two branches represented by equations ( 46 ) and ( 47 ) were 
used. 


To obtain a flow which is symmetric about the line 0=0 as well 
as about the straight streamline the stream function given by 
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$'(0) = ||t(e) + ^(-0^ 


(h8) 


shall he considered. In this formula ^(6) is given hy eqmtion (^7) 
and hy equation (^) if 0 > 0 and 0^ > ^ (-s)3. in the region 


2 , s3/2 

between the characteristic 0 = (-s) 


, ilrC©) and i{r(-0l are hoth 

given hy equation (^7). The choice of the factor 2/3 was made so 
that, in the region of the upper quarter plane 0 > 0 and s > 0 and 

2 , _,3/2 


the section in the left quarter plane where 0 > ^ (-s) ■ 
function is simply given hy 


the stream 


$ = — 




(i*9) 


The physical plane for this portion of the hodograph has already heen 
computed hy Ehlers (reference 13). 


The relation for the stream function in the region between the 
2 
J 


characteristics 0 = is 


where 


^ k /2\1/3, ^i/2\~ 6(0) 6(-0)*l 

|C0S — ( 50 ) 


6 = cos 


-1 0 

I (-s)3/^ 


6(0) = It - 6(-0) 


0 >0 
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Hence, with the aid of a simple trigonometric identity the stream 
function 

‘ 51 ) 

is obtained. The streamlines in the hodograph plane for the entire 
flow given hy equation (W) are shown in figure 11(a). 

To compute the coordinates for the curves of constant speed near 
the y-axls the relation 


z 


P0 


1 

1 


ue=o 



d0 


was employed. By a MacLaurin expansion, for small values of 0, the 
X coordinate can he computed from 


X 





0=0 


$00 + o{03) 


where $q is the value of the stream function when 0=0. The y 
coordinate is computed from 



For small values of 0, this is approximated hy 




where cp = 0 at 0=0. The flow in the physical plane is shown in 
figure U(h). 
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The following paragraphs will give a discussion of the flow pattern. 
The derivatives and fg of the stream function in equation (51) 

2 3/2 

become infinite at the line of branch points 0 ±-^(-s) =0. These 

curves are the characteristics in the holograph plane which eire cusped 
at the origin of coordinates. As these characteristics are approached 
along the branch of the stream function given by equation (49) 
and $g remain finite; thus there is a jump in the first derivatives 

at the characteristics 6 (-s)^^^ = 0. Furthermore along the 

$ = 0 streamline, fg and have the following values: 


^0 


(2/3) 1/3 

o^s 


f S = 0 


The quantity is seen to be singular at the sonic point. 

The effect of the above singularities on the character of the flow 
will now be investigated. From' the equation transforming the solution 
from the holograph to the physical plane, namely. 


(■^"Vccnstant = ^ ^ '1= 

and the relation 

2 

ds = ds* = — dq 

the velocity gradient along a streamline can be expressed in the 
following form: 


^ / ^0 \ 
\di/$=k + ^$0 / 
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where J is the Jacohian of the transformation from the hodograph 
plane to the physical plane. As the Mach line corresponding to the 


characteristic 


e = I (-s)3/2 


is approached from the sonic point of 


the streamline, the velocity gradient approaches a finite value since 
and J are hoth finite. As the same Mach line is approached from 
the direction of the 0 = 0 line, the quantities J — oo and 


^0 - 
from 


00 in such a way. that 


f— ^ 




is finite hut differs in value 


when approached from the "other direction. The line 


J = 00 then corresponds to a discontinuity in the velocity gradient 
along the streamline. For the center streamline f = 0, the velocity 
gradient becomes 


2 


(52) 


Now for $ = 0, $0 ^ as s — * 0. Since sjg has opposite signs 

on the "two sides of the slit corresponding to f = 0 in the hodograph 
plane, and approaches a finite nonzero value as s — > 0, there is a 
jump in the velocity gradient given by 


(: 




(2q*y 


(iA532)V^|J 


1/3 


where q-’^ and p* are the critical velocity and density, respectively. 
To study the curvatiire of the streamlines, the relation 


\l®/^=Constajit 



is considered. In a manner similar to the computations of the velocity 
gradient along the streamline, the following result is obtained for the 
curvature: 



HACA TW 25^7 


39 


This expression is discontinuous at the line J = «>^ since ^ approaches 

different values as the characteristic (J = ») is approached from the 
two directions along the streamline. Thus, the Mach line (J = oo) is a 
curve at which the cuivature of the streamlines as well as the velocity 
gradient is discontinuotis. 


SOME OTHER SOLUTIONS OF THE TRICOMI EQUATION 


Trlcoml (reference 11) has shown that a single-parameter family of 
solutions of the differentleil equation 






is given hy 




(55) 


where a = s3 + 0 “^, Now it is possible to derive from this formula 

some solutions in closed form. In place of the hypergeometric function 
above, consider the following related functions (see p. 286, Whitaker 
and Watson (reference 19)): 


ijr'n = o^( 2o_\“^+V3 j,/ ^ 1 ^ 2 + it. — 2ctJ\ 

\ a + ej V 3 3 3 + 0/ 


and 


t’2,n = * 5= 7f?) <57) 


2 1 

In ilr' let n + — = mj and in ilr’ let -n + — = m. Then the 
^ l,n 3 ^ 2,n 3 


solutions become 
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ho 


x,m 


K” ■ ? ^) 


(58) 


^’2,m = 


7T?) "■(”' - 5' 


Thus it is apparent that 


^ l,m ^ 2,m' 


(59) 


for a given value of m. 

In order to see how these solutions may he summed, consider the 
hypergeometric function written in series form; namely, 


r + m - + n)z^ 


r(m)r 

(m-l) 

r( 2 m) 

CO 

r- 

r(m) 

z_ 

n=0 

r(2m) 

\ 

d 

r(m) 

dz^l , 

r(aa) 

d 


r(2m + n)nl 


n=0 


r + m - + n - l)(m + n - 2) . . . (n + l) 

r ^ - 2)^^ + n - 1)1 


r(m) dz^^-l 


^ni 


C 


Pz 


UO UO, 

— » ' n=0 

2m-l 


r(n + m - i) _ 

\ 1-; ^ 4 d*a.-i 

Z_ 
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The Integrand can be recognized now as simply the binomial expansion of 


/ 


(1 


z) 


l/3-m 


Thus the. expression for and ijf'g becomes 



r 


2,m 


^-2nn-l. ^ 


l^m' 


(61) 


The integrations involved in the above formula can be very easily per- 
formed. Putting m = 1 in equation (60) yields 
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Putting m = 1 into equation (6l) yields 

Putting m = 2 into equation (6o) yields 

^1^2 = - (e - cr)^/^ - cr|Ie + a)^/^ _ (q _ 

Now (0 + and (0 - are separate solutions of the dif- 
ferential equation. The functions ^ 2,1 constructed from 

solutions of the form differentiating with respect to 0. 

Such a process of successive differentiation yields the following solu- 
tions of the Tricomi equation: 


ilr’i = (a + 0 ) 


1/3 


f 2 = 


(g. 0)^/3 


(a +0)^/3 3(ff+0)V3 

V 3 5 i;:' 


o2 g3 


27 (ff + e)'^^3 9 (g + 0)^/3 _ (cr + e)^/^ 

^ ^ = 28 „3 


V ( 62 ) 


cr g 

§i(<7+6)^°/3 . 

5o g' 
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These five fimdamental solutions, together with the set with -0 
replacing 0 may he combined to form the solutions corresponding to 
m = 1 in eqiiatlon (60) and m = 1, 2, 3^ and 4 in equation (6l), 

respectively. By equation (59) > it is seen that cn}r'2, cr^T|r ' 3 j k> 

and cr'T^iJr*^ are also solutions of the differential equation. This may 

he verified hy substituting them into the differential equation. 


The functions and have the property that the 

point s = 0 and 0=0 is mapped into the point at » in the 
physical plane. In fact the solutions all become indeterminate at that 
point. Thus, these solirtions may be used to construct flows having a 
uniform sonic stream at infinity. If 0 is replaced by 0 - iX and 


the real part Is considered, then the point 


0 = 0 and # s 


3/2 , 

' = \ maps 


into infinity in the physical plane. Tomotlka and Tamada used the 
stream function 


'if 


Re 



0 - iX) - ti(-0 + iX) 





to construct the flow over a cusped body having locally supersonic 
regions . 


The solutions o3jt3(-e,s) - +3(6,62!, 

S,b) - ti,.(0,s)j, and cr'^|t5(-0,B) - ^5(6^8^ may be used to 
construct channel flows having locally supersonic regions. Let 


n = jj (0 - ±xf + I s3 - (e - ix^ 

Since 

Ws3n 

J(0 - l\f + ^ s3 + (0 - ix^ ^ 

the appropriate solutions may be written 
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where a(\) = ^(0 - IX)^ + ^ s^. Expressing the real parts in terms of 
the magnitudes and arguments of the complex quantities yields 
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If 6 goes to zero through positive or negative values for 
s > then oi equals 2it or zero, respectively. Hence 


Pi = and the stream functions vanish identically. For 

/o \2/3 ^ fUL 5 7^5 V-/3 

s < Ij , however, cu = it at 0 = 0 and 3i = 


X9 + X 


and the stream functions reduce to non-trivlal functions of s. Thus, 
the stream functions in equations (63)> (64), and (60) will yield 
channel flows with locally supersonic regions. It should he pointed 
out that the solutions given by equation (60) for all integral values 
of m > 1 seem to he reducible to the form 


V = . 0)1/5] 

The solutions for higher values of m than those given in equation (62) 
have a similar form and may he combined to give channels having locally 
supersonic regions. It is then apparent that it is possible to approxi- 
mate quite closely the flow through a given nozzle by taking a sum of 
these solutions, providing the Jacobian of the transformation to the 
physical pleine does not vanish. 

The stream functions given in eqtiations (63)^ (64), and (65) with 
X = 0 cannot be continued into the supersonic region, since, at 0=0 
and s = 0, '^q and ilfg are both zero and hence the Jacobian vanishes 

at this point. These solutions, however, may be combined with the 
solution in equation (l4). 


Brown University 

Providence, R. I., March 13^ 1950 
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TABLE I.- TABULATION OF s, l/q, q, AND X/p 


s 

i/q 

q 

i/p 

0.3^ 

1.45156 

0.68891 

1.27356 

. .26 

1.35002 

.74073 

1.33361 

.18 

1.26118 

.79291 

1. 39801 

.16 

1.24065 

.80603 

1.41545 

.14 

1.22071 

.81920 

1.43333 

.12 

1.20135 

.83240 

1.45175 

.10 

1.18252 - 

.84565 

1. 47078 

.08 

1.16420 

.85896 

1. 49050 

.06 

1.14636 . 

.87233 ■ 

1.51096 

.04 

1.12897 

.88576 

1.53223 

.02 

1.11200 

.89928 

1.55437 

.00 

1.09544 

.91288 

• 1.57744 

-.02 

, 1.07927 

.92655 

1. 60151 

-.04 

1.06345 

.94034 

1. 62666 

-.06 

1.04797 

.95423 

1. 65293 

-.08 

1.03282 

• .96822 

1.68o42 

-.10 

1.01796 

.98236 

1.70918 

-.12 

1.00339 

.99662 

1.73929 

-.14 

.98908 

1. ono 4 

1.77084 

-.16 

.97502 

1.02562 

1.80387 

-.18 

.96120 

1.04037 

I.838W 

-.22 

.93420 

1. 07044 

1.91318 

-.27 

.90152 

1. 10924 

2.01693 

-.32 

.86986 

1. 14961 

2.13391 
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TABLE H.- TABULATION OF s, s*, AND a>]_ 


s 

s* 

0 ^ 

0 . 3 ^ 

0.20075 

1.16534 

.26 

.14508 

1.23307 

.18 

.09521 

1. 30080 

.16 

-08355 

1.31774 

.15 

.07782 

1. 32620 


.07217 

1. 33467 

.12 

.06109 

1.35160 

.10 

.05028 

1.36854 

.08 

.03973 

1.38547 

.06 

.02944 

1.40240 

.oij- 

.01939 

1. 41933 

.02 

.00958 

1. 43627 

.00 

.0000 

1.45320 

-.02 

-.00936 

1. 47013 

-.04 

-.01851 

1. 1*8707 ■ 

-.06 

-.02745 

1 . 504 oo 

-.08 

-.03619 

1.52093 

-.10 

04474 

1.53786 

-.12 

-.05311 

1.551*80 

-.14 

-.06129 

1.57173 

-.16 

-.06930 

1.58866 

-.18 

-.07714 

1.60560 

-.22 

-.09233 

1.63946 

-.27 

-.11046 ■ 

1.68179 

-.32 

-.12077 

1. 72412 
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Flgjure 6.- Mach, number as a function of s*. 
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(a) StreamllneB and constant-velocity lines. 

Figure 8.- Constant-velocity lines in pliysical plane for 
Sq = 0.l8(X = 0.0509117). 
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(b) Physical plane. 


Figure 11.- Concluded 











